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MECHANICS OF AEOLIAN PROCESSES- 
SOIL EROSION AND DUST PRODUCTION 

by 

M. M. Mehrabadi 

Dept, of Mechanical Engineering 
Tulane University 
New Orleans, Louisiana 70118 


ABSTRACT 

Aeolian (wind) processes occur as a result of 
atmosphere / land-surface system interactions. A 
thorough understanding of these processes and their 
physical/mechanical characterization on a global 
scale is essential to monitoring global change and, 
hence, is imperative to the fundamental goal of the 
Earth Observing System (Eos) program. 

Soil erosion and dust production by wind are of 
consequence mainly in arid and semi arid regions 
which cover 36% of the Earth's land surface. Some 
recent models of dust production due to wind erosion 
of agricultural soils and the mechanics of wind ero- 
sion in deserts are reviewed and the difficulties of 
modeling the aeolian transport are discussed in this 
report . 
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MECHANICS OF AEOLIAN PROCESSES - 
SOIL EROSION AND DUST PRODUCTION 

1 . INTRODUCTION 

Aeolian (wind) processes occur as a result of 
atmosphere/ land-surface system interaction. Study of this 
interaction is one of the broad geologic science issues 
identified by the Earth Observing System (Eos) Science and 
Mission Requirement Group (Butler et al., 1984) and the 
Eos Science Steering Committee (Butler et al., 1987). A 
fundamental goal of Eos is to develop an improved under- 
standing of the processes that control the formation and 
evolution of the solid Earth. Aeolian processes are per- 
sistent sculptors of the face of the planet Earth (and 
Mars) . Hence a thorough understanding of these processes 
and their physical/chemical characterization is essential 
to monitoring global change and, thus, is imperative to 
this goal. 

Wind erosion and transport are active geologic 
processes that occur mainly in arid and semi-arid regions 
which cover 36% of the Earth's land surface. Arid and 
semi-arid areas suffering from "desertification" undergo 
destruction of native vegetation which is quickly followed 
by an excessive soil erosion. Examples of areas with 
severe soil erosion are northern China, including part of 
the Loess Plateau and Inner Mongolia, the Thar desert of 
eastern Pakistan and northwest India, and the U.S. Great 
Plains, including parts of Texas, Oklahoma, and New 
Mexico. 

Although soil erosion by wind is of consequence 
mainly in arid and semi-arid regions, it may occur 
wherever soil, vegetation, and climatic conditions are 
conducive. Generally, wind erosion is effective when soil 
is loose, fine, and dry, and has little or no vegetation. 
The present rate of fertile soil loss in the United States 
may be as high as 100 tons per acre per year in half a 
dozen major agricultural regions, including the corn-belt 
states of Iowa, Missouri and Illinois. The Council on En- 
vironmental Quality (Sheridan, 1981) estimates that about 
10 percent of the U.S. land mass is in a state of severe 
or very severe desertification. 

Large dust storms originate in major deserts of the 
world (Fig. 1). In most of these deserts, dust is formed 
by in situ weathering of bedrock and from abrasion of the 
bedrock by the wind-blown sand (Goudie, 1978). Large, 
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frequent dust storms originate from normally semi-arid 
areas (e.g., the Great Plains of the United States; 
Central India; and the Russian steppes) that periodically 
become arid, undergo abnormally strong windy periods, or 
have their vegetation removed by man or nature (Pewe, 
1981) . 

In terms of particle size, dust can be divided into 
two major groups depending on how far a distance they can 
be transported by the wind (Pewe, 1981) . Dust particles 
that are carried between a few kilometers to less than 100 
km, by dust devils and dust storms, are mostly between 
0.005 mm to 0.05 mm. The other type of dust consists of 
particles between 0.002 mm and 0.01 mm in size. This 
fine-grained, sorted dust is carried to the troposphere as 
an aerosol (i.e., a dust that remains suspended in the air 
until brought down by rainfall onto water and land 
surfaces) . Dust from the Sahara collected in the Caribbean 
(Fig. 2) has been reported to be more than 98% smaller than 
0.01 mm (Prospero, et al., 1970). 

Aerosols affect the reflectance, transmittance and 
absorptance of the atmosphere (Verstraete, et al., 1988). 
They modify the amount of shortwave radiation absorbed and 
longwave radiation emitted by the climatic system. Their 
impact on the climate is further enhanced by their role as 
cloud condensation nuclei. 

Aerosols impact on the transfer of radiation through 
the atmosphere also affects the interpretation of the 
satellite remote sensing data. While this creates a 
nuisance for data obtained using passive remote sensing 
techniques, aerosols provide excellent tracers when active 
instruments such as the Laser Atmospheric Wind Sounder 
(LAWS) are used for wind velocity measurements. 

Dust has many other serious environmental conse- 
quences. These include its effect on crop growth, ocean 
sedimentation and its effect on the acid/base balance of 
atmospheric deposition. Dust storms cause disease, suf- 
focation of cattle, development of static electricity, 
disruption of transportation, and destruction of property. 

A model for the estimation of total dust production 
for the United States has been discussed recently by Gil- 
lette (1988). After reviewing this model in Sec. 3, the 
mechanics of aeolian transport of particulate materials 
is discussed in Sec. 4. The long and short term goals of 
this research effort are described in the next section. 


2. OBJECTIVES 

Some recent literature on mechanical modeling of 
aeolian processes are reviewed in this report. Specifi 

cally, two kinds of models are discussed. The n^nroduc- 
deals with the agricultural soil erosion and dust produc 
tion by wind. The second kind deals with the more bas c 
problem of wind erosion in deserts. 

This effort is part of a comprehensive program of 
research focused on assessing the potential impact of 
desertification on climate changes and global 

habitability. The ultimate goal of this ^mechanicl of 
develop a quantitative understanding of the ° f 

atmospheric coupling with land ; systems as rela ^® d 

desert processes in and and semi-arid regions. Exist g 
illl banks , data generated by laboratory experimentation 
and, later, data generated by Eos sensors will be used in 
the validation stage of the model. 

The short term goal of this research is the formula- 
tion and validation of a model for the main mechanisms in 
volved in the aeolian transport of particles based on 
principles of mechanics at the particle-scale. A 

present, no comprehensive quantitative analysis of such 
processes is available. 
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3. A MODEL FOR DUST EMISSION BY WIND EROSION 

A model for the total dust production for the United 
States, recently proposed by Gillette and Passi (1988) is 
reviewed in this section. The primary use of the model is 
m the inventory of alkaline elements for use in acid/base 
balance studies of atmospheric precipitation by the Na- 
tional Acid Precipitation Assessment Program (NAPAP) . 


Wind Profile in the Lower Planetary Boundary Layer 

In the atmospheric inertial sublayer, which is the 
lower part of the planetary boundary layer, the wind speed 
can be assumed to vary with height according to the 
logarithmic wind profile (see, e.g., Priestely, 1959, or 
Plate, 1971) 


U(x} 




( 3 . 1 ) 


if the stratification is nearly adiabatic ( for winds 
above the velocity required for wind erosion of soils, 
Gillette, 1981, reported that the atmosphere was neutrally 
stratified). In (3.1), z # denotes atmospheric roughness 
height, taken to be 1 cm for smooth surfaces (e.g., air- 
port surfaces) , and V* is the friction velocity defined 
m terms of the shearing stress (shearing force per unit 
area of ground surface parallel to the wind direction) as 



or 


tmy ; 1 




where $ is the air density. Note that V* is the velocity 
of the wind at a height approximately equal to 1.492 z 0 . 

It is convenient sometimes to introduce a drag coef- 
ficient, Cj , defined by (Priestely, 1959) 



where (3.3) 4 has been written using (3.2) x to eliminate x. . 
Gillette and Passi (1988) use a method proposed by 
Wieringa (1976) to evaluate the drag coefficient Cj . 
This procedure is summarized in the next two paragraphs. 
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Wieringa (1976) takes the wind speed at a height of z w 
= 60 m to be regionally constant. In other words, he as- 
sumes that the wind speed at that height, U(z^), does not 
show excessively the influence of the particular roughness 
element close to the wind station where an average wind 
speed U ( z 3 ) is measured at a height z, . Suppose that at 
the same location we wish to know the average wind speed 
U(z r ) over a hypothetical flat open reference terrain at a 
height z r . Referring to Fig. 3, and calculating the 

slopes of the two wind profiles, will give 

j* L* a Vz»» fj.V) 

uCzO-UCZs) * 

and 


Uiz k )-U(2 r ) 


Z 7Zor 

u( Zr) 


Eliminating U(z fe ) between (3.4) and (3.5), one finds 


UUr X ; - Ufz.'l 


u Vx. 


M 


This is the required relation between U(z r ) and U(z s ). 

If z r is taken to be 1.492 z.r , so that In z r /z 6r = 
0.4, then from (3.1) it follows that U(z r ) = ^ and from 
(3.6) it follows that 


V, - U(z s ) -11 — 




** Vz., 2 Vz #r 


Using (3.7), the drag coefficient can be calculated from 
(3.3) as 

_ «<■ 1 * 

c ‘ * W. Ty M 

I 'Znr 
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Gillette and Passi (1988), use a value for z 0 , - 1 cm, 
typical of clipped grass at airports, and z^ - 0.002 cm, 
from an average measured aerodynamic roughness heights 
they obtained in many experiments. With these values they 
find that 



( 3 . 1 ) 


where z s is in cm. Substitution of (3.9) in (3.3) yields 

V. --£i 2- U rt.iol 

which relates the friction velocity to the wind speed. 


3.2 Threshold Friction Velocity 

Threshold friction velocity for wind erosion cor- 
responds to the minimum wind stress needed to overcome 
forces holding soil particles in place. Threshold fric- 
tion velocity is denoted by V*t and is defined in terms of 
surface shear stress at threshold (or "threshold wind 
stress"), , by the relations 



Although one must rely heavily on the values of the 
threshold friction velocity for estimation of dust produc- 
tion and wind erosion, there exists no theory for relating 
this quantity to the physical characteristics of the soil. 

Threshold friction velocity has been determined ex- 
perimentally for idealized surfaces by many investigators 
(see Gillette, 1988, for references) starting with the 
pioneering work of Bagnold (1941) who proposed a formula 
for V*t in terms of the particle diameter and density (see 
section 4.1 of this report). The physical circumstances 
in agricultural soils is complicated, however, because of 
a wide variety of particle sizes and the effects of 
w ®tting/ drying , freezing/thawing, vegetation, and aggrega- 
tion of particles. Measurements of the threshold friction 
velocities for dust production for a wide class of 
agricultural soils has been reported only recently by Gil- 
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lette (1988) . Gillette and Passi's estimation of dust 
production is based on the measurements of Gillette (1988) 
and the values of threshold friction velocity obtained for 
arid regions by Gillette et al. (1980, 1982) . 


3 . 3 Dust Flux as a Function of Wind Speed 

-> 

Figure 4 shows the vertical flux of dust, G (g cm s 
) as a function friction velocity. The data were reported 
by Gillette (1974,1981) who measured dust fluxes in out- 
door eroding fields, Fairchild and Tillery (1982) and Bor- 
rman and Jaenicke (1987), who measured fluxes with indoor 
wind tunnels. 

Gillette and Passi (1988) express the flux of dust, 
G, by the relation 


G 



CJ./JL) 


where C, is a constant. This is shown as curve A in the 
Figure. Gillette and Passi (1988), state that (3.12) is 
based on an unpublished theory of Owen (personal com- 
munication between Gillette and Owen, 1987) . 

Substituting from (3.10) into (3.12), gives 


s(u-) - 



( 3 ' 3 ) 


where U t is the threshold wind velocity obtained from 
(3 . 10) , i.e. , 



0.13 

U 2 S 


«L 


(i-H) 


Let p (U) be the probability density function of the 
wind speed during the time period of interest. A widely 
used probability density function for natural winds is the 
two-parameter Weibull distribution (Gillette and Passi, 
1988) 


P(U) 



( 3 . 15 ) 
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where the two parameters c and k are determined directly 
from the wind data. The data set used by Gillette and 
Passi (1988), "the Wind Energy Resource Information System 
(WERIS) " , provided monthly values of the mean hourly wind 
speed and the "pattern factor" statistics 

X =» frtt) 

II* 

where a bar over a quantity indicates the mean of the 
quantity for each month. WERIS provided data for a 31- 
year period, 1948-1978, for 1432 wind measuring locations 
in the United States. It turns out that the pattern fac- 
tor x is a unique function of Weibull shape parameter, k, 
so that 


rp-OA) 

r’o+'V) 


(3J7) 


where T is the gamma function. Given x, (3.17) can be 
solved for k by trial and error. 

With the probability density - function of the wind 
speed p(U) , given by (3.15), G(U) is integrated for all 
wind speeds above the threshold, CL for each location 
(subarea) of the region of interest, l.e., 

/ G(u) f (u)JU 

U * (*«) 



where G(U t ) is the average value of the flux estimates, 
G (U) , for one sub-area, and where (3.13) and (3.15) are 
used to write (3.18) 
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3.4 Estimation of Total Dust 

The region of interest is divided into N sub-areas. 
Then the average value of the flux estimate for each sub- 
area G(U t ) (mass per unit area per unit time) is calcu- 
lated from (3.18). This averaged estimates are then 
summed over all erosion areas within the region of inter- 
est, using various parameters as weighting function. Thus 
the mass of dust, E, emitted in the time period &T is 
written as (Gillette and Passi, 1988) 

N 

£=C 21 $(Lj) A; G ( U t; ) AT 

where C is a constant to be determined by calibration (the 
constant C, coming from G is lumped into C) ; i is the index 
of summation over N different erodible areas within the 
region of interest; R ; is the effect of soil roughness; 
g(Li) is the effect of field length, is the area of 
the land being considered. 

Roughness of the soil has the effect of trapping 
soil particles, increasing the friction velocity of the 
wind and thus inhibiting erosion. For this reason the ex- 
act effect of soil roughness is complicated. Gillette and 
passi (1988) examined data from the U.S. Department of 
Agriculture 1982 National Resource Inventory (NRI) for 
soil information, land use, and wind erosion parameters to 
evaluate a "ridge roughness function" proposed by Armbrust 
et al. (1964) which was based on wind tunnel tests. Since 
over 97% of these values examined in the data set were be- 
tween 0.5 and 1 (for prominent erosion areas located 
largely in the panhandles of Texas and Oklahoma) , they 
conclude that the maximum effect on dust production can 
not be more than a factor of 2. 

The effect of field length on dust production is 
thought to be (Gillette and Passi, 1988) (i) an 
aerodynamic feedback mechanism whereby the airborne sand 
increase the aerodynamic roughness height by absorbing 
forward momentum from the wind and thus decreasing its 
speed, and (ii) an increase of soil flux with increasing 
field length due to sandblasting of the loosened downwind 
soils. Gillette and Passi (1988) assume g(L;) to be 1, 
ignoring the effect of change with distance downwind of 
the dust flux. It should be noted that mechanisms 
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similar to those described in (i) and (ii) are incor- 
porated in a recent numerical model considered by Anderson 
and Haff (1988) . 

To evaluate C in equation (3.19), Gillette and Passi 
(1988) make use of data from total dust production for a 
large area including the panhandles of Texas and Oklahoma. 
They find a value of 1.4 x 10 _w for the constant C. 
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4. AEOLIAN TRANSPORT OF SAND IN DESERTS 

A major contribution to the study of the physics of 
wind-blown sand in deserts was made by the classic work of 
Bagnold (1941) who combined extensive in-situ observations 
and wind-tunnel experiments to describe the motion of sand 
particles in the air. 

Transport of dust, sand and other particulate 
materials by wind takes place mainly in the form of sur- 
face creep, saltation, and suspension (Bagnold, 1941) . 
Surface creep, also known as ’'traction" is the movement of 
particles along the surface. Saltation refers to a leap- 
ing motion of particles where they lift off the surface 
and travel in ballistic-type trajectories in the air 
before returning to the surface. Suspension is the motion 
of (less massive) particles that are carried upwards by 
the wind to travel at indefinite heights within atmos- 
phere ( Fig. 5). 

Bagnold (1941) states that of the total sand in mo- 
tion in air about three-fourths move in saltation and a 
quarter in surface creep. The surface creep consists of 
slow jerky advance of the surface grains which are knocked 
along the surface by the impact of descending saltating 
grains. According to Bagnold, a saltation of fine grains 
can maintain a surface creep over a bed composed of grains 

far too large to be moved by the direct action of the 
wind. 


Bagnold notes that a flat sand surface is inherently 
unstable. He argues that since the grains are of varying 
size, the surface is not perfectly even and tiny depres- 
sions can form such as the one shown, much magnified, in 
Fig. 6. The series of parallel lines shown in the figure 
are supposed to represent the descending saltating grains. 
On the lee side, along the line AB, the points of impact 
are far apart, but on the windward side many impacts oc- 
cur. This means that the original depression will get 
bigger. Eventually, grains accumulate at C and a second 
lee slope CD is formed, and so on. Using this line of 
reasoning, Bagnold explains the formation of ripples and 
argues that the range of the characteristic path of grains 
is coincident with the ripple wavelength (Fig. 7) . 

Since saltation is the main mechanism of aeolian 
transport in sand, the remainder of this section is 
devoted to a review of some recent works on this topic 
(White, 1986, Anderson and Hallet, 1986). 
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4 . 1 Threshold Conditions 


Bagnold (1941) has shown that at the initiation of 
saltation, the surface shear stress at threshold is 

% - co.o 


where A is a varying empirical coefficient, is the par- 
ticle density, g is the acceleration due to gravity, and D p 
is the mean particle diameter. In terms of the friction 
speed at threshold, V #t , 



ft.4) 


Since the particle density (2.64 g/cm* for quartz) is much 
greater than the atmospheric air density (1.22 x 10"* ) we 
have 


if. 

S 


He Altv » I 


Thus (4.2) can be written as 





A 


*y%P? 




Bagnold and others have shown that A = A(B), where 




( 4 . 5 ) 


is the particle friction Reynolds number. In (4.5), * is 
the kinetic viscosity of air. The variation of A with B 
obtained from wind tunnel data is shown in Fig. 8 (after 
White, 1979). Note that A « 0.118 for values of B > 10. 
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4.2 Forces acting on the Particle 

Forces acting on the particle are (i) body force 
which resnits from gravity, F & , (ii) aerodynamic drag, f„ 
' aerodynamic lift, f l , and (iv) Magnus lift which° 

results from particle spin, F M . The trajectory of a typi- 
cal saltatmg is schematically represented in Fig. 9 . Note 
that a number of forces including the foce due to buoyancy 
( since 3 << ) and the interparticle force have been 

neglected; see White (1986) for a detailed discussion of 
these and other forces. 

. , . , ro dynamic drag is a function of the velocity of the 

particle relative to the air. v 


V . 


U-V 
* Zp 


(♦.*) 


where (with components x and z ) is the particle 

velocity. The drag force is expressed as 




(V.T) 


where C b is the drag coefficient which is a function of 
the particle Reynolds number. Re = v^Dp/v, i.e.. 


C * - C t>( *0 


C*.i) 


lift is due to shear flow and its mag- 
nitude is given by (see Anderson and Hallet, 1986) 


w c l(^)»(«C-«C) 


C^-4) 


and U.*, are the air speeds, at heights cor- 
responding to the top and bottom of the grain, and C, is 
given by (Chepil, 1958) y ' •- ls 


'L * °‘SSC D 




The Magnus lift is given by (see White, 1986, and An- 
derson and Hallet, 1986) ' ana An 


XX-13 



~h 

where A is the angular velocity of the particle. 

Rubinow and Keller (1961) derived the following rela- 
tions for the moment, M, on a particle due to its spin XL , 

where f is the air viscosity. Given the forces and the mo- 
ment acting on the particle, equations of motion for the 
particle can be written and a numerical solution for the 
particle trajectories can be obtained (see White, 1986, 
and Anderson and Hallet, 1986) . 


4.3 Saltation Trajectories 

Given a liftoff angle, ot , speed of liftoff Upo , and 
initial particle spin, -fl.© , the particle trajectory is 
numerically obtained using time iterations. At each time 
step, the Reynolds number is calculated, a drag coeffi- 
cient is found, and the spin rate and local shear are cal- 
culated. Using these, accelerations, velocities, and dis- 
placements are updated for each increment. Results ob- 
tained by White (1986) are shown in Fig. 10. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


Some recent models of dust production due to wind 
erosion of agricultural soils (Gillette and Passi, 1988) 
and mechanics of wind erosion in deserts (Bagnold, 1941, 
and others) were reviewed in sections 3 and 4 of this 
report. It can be seen that in spite of the progress made 
so far, there are several issues that need more attention. 

within the framework of the models described two 
specific issues should be addressed as a next step (i) 
modeling of the threshold velocity particularly for 
agricultural soils, relating it to texture, moisture, 
ere,, and (ii) the initiation of saltation which is still 
a source of controversy. The actual mechanism of particle 
liftoff and the forces operating at the initial stage are 
not well understood. The recent work of Anderson and Haff 
(1988) seem to be promising in this regard. 

More generally, however, there are several viewpoints 
for modeling of particle-gas flows, namely, the "discrete 
particle" or "tracking" approach (see, for example, Durst, 
et al., 1981) and the "non-discrete" or continuum approach 
( see Ahmadi, 1986). The suitability of these approaches 
for modeling aeolian transport of soils must be examined 
in detail. The most significant aspect of the problem is 
modeling the interaction between the air-particle mixture 
and the boundaries of the flow. There is a continuous ex- 
change of mass, momentum and energy between the bed and 
the air-particle system. Unless we take these exchanges 
into account, using basic principles, we have to rely on 
empirical formulas for dealing with wind erosion and dust 
production problems. 
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